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Abstract. The modern study of the dynamics of stellar systems requires the use 
of high-performance computers. Indeed, an accurate modelization of the structure 
and evolution of self-gravitating systems like planetary systems, open clusters, 
globular clusters and galaxies imply the evaluation of body-body interaction over 
the whole size of the structure, a task that is computationally very expensive, in 
particular when it is performed over long intervals of time. In this report we give a 
concise overview of the main problems of stellar systems simulations and present 
some exciting results we obtained about the interaction of globular clusters with 
the parent galaxy. 
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1. Introduction 

The question "What is a supercomputer?" 
has a surprisingly trivial answer in the 
CRAY web site: "A supercomputer is de- 
fined simply as the most powerful class of 
computers at any point in time." Actually, 
what modernly distinguishes a supercom- 
puter from normal computers is its more 
or less complex architecture that deter- 
mines its effective power. Usually, the per- 
formances of the high-class processors are 
quite similar and so producers, in order 
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to increase significantly computing speed, 
assemble an enormous number of proces- 
sors linked in the most efficient way. In 
many cases a relatively cheap solution is 
that called "Beowulf, i.e. the collection 
of a certain number of, say, commercial 
Intel processors linked each other with fast 
Ethernet. A standard type of "Beowulf is 
the so-called Linux cluster. These systems, 
are well suited to a wide range of appli- 
cations that not require top-performances 
but, rather, for them suffices a machine 
whose use is reliable and, possibly, dedi- 
cated. Whoever is experienced in numer- 
ical simulations knows very well that the 
largest quantity of time is spent in testing 
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the code in conditions as similar as possible 
to the ones of the real " long lasting" scien- 
tific runs. It is, thus, evident the advantage 
of having the full compatibility between the 
testing machine and the running machine 
that is guaranteed by using the same ma- 
chine. Moreover, in this way all the prob- 
lems of large data file remote downloading 
is avoided. 

In spite of these considerations, the "per- 
sonal" Beowulf solution is far from be- 
ing exhaustive. Many scientific problems (I 
would say the highest-level ones, i.e .the 
most interesting) require, to be approached 
with an acceptable level of approximation, 
higher performances than those reached 
by "Beowulfs", both in terms of comput- 
ing speed and memory (Linux clusters, 
also when composed by an assembly of 
bi-processor machines, are limited in the 
shared memory storage). For this class of 
problems, where I pose many astrophysi- 
cal problems, the availability of huge par- 
allel main frames is essential. A clear exam- 
ple of this is given by the recent construc- 
tion of the Japanese "Earth simulator", 
a supercomputer explicitly aimed at the 
detailed simulation of wheather evolution 
over a wide space and time domain with an 
unprecedented resolution. This computer 
(resulting first for the second consecutive 
year in the Top 500 international rank- 
ing list http://www.top500.org) is com- 
posed by 640 nodes of 8 vector proces- 
sors each giving a total peak power of 35.8 
Teraflops/s. This enormous power is largely 
due to the use of vector processors, which 
represents a constructive choice different 
from that, most widely used in U.S., of as- 
sembling as many " commercial" processors 
as possible. 

In any case, whichever the supercom- 
puters, it is ascertained that science of ex- 
cellence in general, and astrophysics and 
cosmology in particular, requires them. . 

2. Supercomputing and astrophysics 

Modern astrophysics and cosmology are 
characterized by dealing with complex 



problems whose dynamical range covers 
extended space-time scales. Some specific 
topics where this happens are (from small 
to larger systems): 

• stability of solar system, 

• star formation, 

• final phases of stellar evolution, 

• dynamics of galaxies and stellar sys- 
tems, 

• physics of AGNs and Quasars (mass 
accretion onto compact objects), 

• large scale structure. 

A correct comprehension of these top- 
ics requires the analysis of the contribu- 
tion of the many processes involved, act- 
ing from the micro space-time scales (e.g. 
atomic and molecular heating and cooling 
processes in star formation or relativistic 
supersonic accretion of matter onto black 
holes in galactic nulcei activity, etc.) up to 
the macro structure of the system mainly 
governed by gravity which is, in addition, 
the main engine of activity also on the 
micro-scales, due to the irreversible trans- 
formation of gravitational energy into ther- 
mal. 

It is nowadays absolutely clear that 
firm conclusions about the mentioned open 
problems cannot be reached by means of 
simplified models. Consequently, the only 
way to follow in detail the mutual interac- 
tion among phenomena acting on different 
space-time scales up to the clear vision of 
the overall evolution of the system relies 
on the production of sophisticated numer- 
ical modelizations. Due to the complexity 
of problems, running these models require 
indeed high performance computers, i.e. su- 
percomputers. 

2.1. The role of gravity in numerical 
astrophysics 

By the point of view of computational 
weight, the most important difference be- 
tween astrophysical and terrestrial fluid- 
dynamics is the presence of self-gravity. 
This implies the substitution in the equa- 
tions of motion of the self-gravity body 
force to the simple fixed force field. 
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Table 1. Stellar systems time scales and other characteristics relevants to the types of 
simulations needed. The age of the system is t. In the last row F-P means Fokker-Planck 
method, tree-c. means tree-based codes, PM and P3M stand for Particle-Particle and 
Particle-Particle-Particle-Mesh algorithms, respectively. 



stellar system: 


open cluster 


glob, clus., gal. nucleus 


galaxy, clus. of gal. 


N: 


< 10, 000 


10 b -5- 10 u 


>10 1U 


gravity: 


newtonian 


newtonian, gen. rel. 


newtonian, gen. rel. 


time scale ordering 


trel ^ icr ^ t 


tcr ^ irel ^ t 


tcr ^ ^ t *C trel 


regime: 


collisional 


secularly collisional 


collisionless 


technique: 


gas+direct N-body 


F-P,dircct N-body,tree-c. 


trec-c.,PM,P3M 



Computationally, this causes an enormous 
overcharge, due to the long range charac- 
ter of gravity. This problem is not just a 
prerogative of collision-dominated systems 
(fluids, where the relaxation time scale is 
much shorter than any other dynamical 
time scale) but typical also of collisionless 
or partially collisional systems, like galaxies 
and/or stellar clusters. The long-range be- 
haviour of gravity implies that one cannot 
neglect the interaction with distant bod- 
ies, simply because their number at a fixed 
distance increases, and this corresponds to 
an obvious heavy computational weigth at 
any time step. Things are even more com- 
plicated because of the divergence of the 
newtonian gravitational potential for mu- 
tual distances approaching to zero. This is a 
serious problem whenever fluctuations over 
the mean field are large in small stellar sys- 
tems; the short-range behaviour of gravity 
may induce a dramatic reduction of time 
stepping when a close encounter between 
two, or more, "particles" makes their accel- 
eration to increase up to values such that an 
acceptable error in the integration of their 
trajectories requires a much finer time ad- 
vancing. 

2.2. From few- to many- body systems 

The study of the stability of solar system 
is a typical few-body classical gravitational 
problem. 

The number of bodies (planets, asteorids, 
etc.) is small and the initial conditions rel- 



atively well known so to avoid the prob- 
lems caused by the double singularity de- 
scribed above. By the way, the secular sta- 
bility of the solar system is still an unan- 
swered problem due to the intrinsic diffi- 
culty in evaluating the role of the reso- 
nances that determine the long time be- 
haviour of the system. A correct treatment 
cannot neglect the planet structure, be- 
cause tides favour resonances among plan- 
ets and satellites. So, powerful computers 
are needed to follow the time evolution of 
a system over very many orbital periods of 
the components with an accuracy such that 
one can trust into results that are heavily 
affected by numerical error accumulation. 
It is intriguing, in this context, that the 3:2 
Neptune-Pluto orbital resonance has been 
found as a result of primitive supe rcomput- 
ing bv lCohen fe Hubbard I (|1965j) . 
The actual check of the valid- 
ity of the resu lts ob tain ed by 

ISussrnan fc Wisdom I fmsft. Ib^kaTl 

(1989) and ISussrnan fc Wisdom I (1992), 
that suggest the solar system being chaotic, 
will probably require the next generation 
of specifically dedicated supercomputers. 

Let us say few words about larger N- 
body problems. 

By "intermediate" N-body systems we in- 
tend stellar systems composed by up to 
few times 10 6 stars, i.e. open and globu- 
lar clusters and clusters of galaxies (in this 
case, of course, the individual body is a 
galaxy) . The study of the evolution of these 
systems is a typical "multiple time scale" 
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problem, where the fine grain of the system 
determines a multiplicity of individual time 
scales (the two-body fly-by characteristic 
times At jj ) that combine with the crossing 
time t cr , which, in its turn, depends on the 
coarse grain structure. The dynamic range 
of an accurate time-stepping would span 
about 5 orders of magnitude. The use of 
time step as the minimum over this range 
would cause a CPU time requirement ex- 
cessively large when the aim is to simulate 
the dynamical evolution of the system for 
a sufficiently long physical time (moreover, 
the error cumulation would be a great prob- 
lem). The way commonly adopted to keep 
accuracy within an acceptable CPU time 
consumption is the use of individual (and 
variable in time) time stepping with accu- 
rate sinchronization. 



Individual time stepping, inserted in a 
clever "tree algorithm" , i.e. such that the 
treatment of small scale fluctuations of the 
force field is done by mean of "local" di- 
rect summation has a great power in deal- 
ing with star cluster dynamics.. Actually, 
both high accuracy and speed are guar- 
anteed making possible reliable simula- 
tions of stellar systems of the size of real 
globular clusters, with a 1:1 representa- 
tion of real stars with simulating parti- 
cles (that means using N up to 10 6 in 
the code). Such performances are reached, 
of course, b y mean of parallel implemen- 
tation s fsee iMiocchi fc Capuzzo-Dolcetta I 
(2002). It may be worth stressing that at 
present, globular clusters are the largest 
systems that can be studied with a 1:1 rep- 
resentation, i.e. with an almost complete 
reliability of the results, contrarily to what 
happens with cosmological and large scale 
simulations, where the ratio between the 
number of real objects (stars, galaxies, clus- 
ters of galaxies) and simulating particles 
(iV) is still very large, making the results, 
that are also affected by finite size effects 
and by the smoothing length of the inter- 
action potential, not yet completely reliable 
and of not easy interpretation. 



3. Globular cluster-galaxy interaction 

While the study of the very long-term evo- 
lution of a popolous globular cluster is not 
yet feasible, because it requires an inte- 
gration over a time of the order of the 
evaporation time, i.e. about one hundred 
times the 2-body relaxation time, that is 
t re i ~ l,500t cr for a 10 6 stars cluster, it 
is nowadays possible to examine with good 
accuracy the evolution of a globular cluster 
in motion in the galactic field. 
Actually, useful information on the mutual 
feedback of the cluster with the galactic en- 
vironment is obtained by simulations which 
cover an interval of time of 10 orbital times. 
We are presently working in the program 
of studying the coupled orbital and inter- 
nal evolution of massive globular clusters 
in the field of an elliptical galaxy, mod- 
eled with a self-consistent triaxial galacti c 
model lCanuzzo-Dolcetta fc Vicari. I lf20031 . 
The main aims of this program are: 1) 
to understand why and how tidal tails 
are formed behind globular cluster mov- 
ing in the galactic field; 2) to evaluate the 
amount of mass loss from clusters to the 
halo-bulge; 3) to check, in a straightfor- 
wa rd way, the hypoth e sis ra ised time ago 
by ICapuzzo-Dolcetta I l|l993ft that a com- 
pact object at the galactic centre may ac- 
crete a significant quantity of mass by mean 
of orbitally decayed globular clusters. 

Here we report just few prelim- 
inary results concerning point 1) 
(they will be extensiy el y pre sented in 
ICanuzzo-Dolcetta et al. I ^2003)) . Figure 
n shows the time evolution of a 3 x 10 s 
solar masses globular cluster made by 
N = 165,000 stars whose masses are 
distributed according to a Salpeter's mass 
function, with lower and upper mass 
cutoffs 0.1 M & and 1.2 M©, respectively. 
The cluster moves on a quasi-circular 
orb it of 500 pc rad i us in the (y,z) plane of 
the ISchwarzschild I l|l979l) triaxial galactic 
model of axial ratios 2:1.25:1. It is well 
evident the development of tidal tails 
along the orbit, that stabilize their pattern 
in the direction of the (clockwise) motion 
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Fig. 1. Left panel: The first period of the clockwise motion of an 3 x 1 5 Mfn globular 
cluste r along a quasi-circular orbit of 500 pc radius on the y-z plane of the ISchwarzschildl 
( 1979) triaxial model (the inner box is a zoom of one of the cluster configuration; in green 
the stars physically belonging to the tails. Right panel: as in the left panel, but referred 
to the second orbital period (snapshots are labeled with time); clumps are evident in the 
arms 



after less than one orbital period (about 
7 Myr). Note that the front arm points 
in the direction of motion, exactly as it 
is seen in the globular cluster Palomar 
5 obtained by lOdenkirchen et al. I l|2002j) 
by mean of wide-field photometric data 
from the Sloan Digital Sky Survey (SDSS) 
(Fig. 0). Another interesting feature that 
has an observational counterpart in the 
same cluster is the formation of clumps 
along the tails (see Fig. fl] right panel). 
Clumps correspond to zones where the 
tails change directions, probably due to a 
non-linear cumulative effect of the tidal 
field. On another side, with regard to the 
inner stellar distribution, we see that inner 
relaxation is speeded up by the external 
field, in the sense that mass segregation is 
enhanced, thus inducing a decrease of the 
velocity dispersion in the inner region in 
ag reement with kinem at ic dat a of Palomar 
5 lOdenkirchen et al. I l|200ll) that show 
that Palomar 5 is a "cold" system. 

All these results are preliminary and de- 
serve a careful confirm and generalization; 



by the way we are convinced that these 
types of numerical simulations are of fun- 
damental importance to give answers to the 
initial evolution of galaxies as due to mu- 
tual feedback among the various primordial 
stellar populations and also to the still un- 
known mechanisms of mass accretion onto 
the AGN engines. 
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Fig. 2. From the SDSS Press release june 3, 2002 Ipdenkirchen et al. I l)2002fi ): color- 
coded map of the distribution of stars emerging from the star cluster Palomar 5 (white 
blob). The two long tidal tails (orange) contain 1.3 times the mass of the cluster and 
delineate its orbit around the Milky Way (yellow line). 
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